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Abstract. A number-conserving cellular automaton is a cellular au- 
tomaton whose states are integers and whose transition function keeps 
the sum of all cells constant throughout its evolution. It can be seen as a 
kind of modelization of the physical conservation laws of mass or energy. 
In this paper, we first propose a necessary condition for triangular and 
hexagonal cellular automata to be number-conserving. The local transi- 
tion function is expressed by the sum of arity two functions which can 
be regarded as 'flows' of numbers. The sufficiency is obtained through 
general results on number-conserving cellular automata. Then, using the 
previous flow functions, we can construct effective number-conserving 
simulations between hexagonal cellular automata and triangular cellular 
automata. 
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1 Introduction 

A number-conserving cellular automaton (NCCA) is a cellular automaton (CA) 
such that all states of the cells are represented by integers and the sum of the 
numbers (states) of all cells of a global configuration is preserved throughout 
the computation. It can be thought as a kind of model of physical phenomena 
as, for example, fluid dynamics and highway traffic fiow [7] and constitutes an 
alternative to differential equations. 

There is a huge literature published in the domain which witnesses the great 
interest in number-conserving cellular automata which gathers together physi- 
cians, computer scientists and mathematicians. Actually, this particular model 
of CA applies to phenomena governed by conservation laws of mass or energy. 

Boccara et al. |1 studied number conservation of one-dimensional CAs on 
circular configurations. Durand et al. |2I3| considered the two-dimensional case 
and the relations between several boundary conditions. These results are very 
useful for deciding whether a given CA is number-conserving but do not help 
much for the design of NCCAs with complex transition rules. 

As for the rectangular von Neumann neighborhood case |5l8j , several nec- 
essary and sufficient conditions to be number-conserving are shown. According 
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to these conditions, the local function of a rotation-symmetric NCCA is ex- 
pressed by the sum of arity two functions as in |4]. Designing the functions, we 
constructed several NCCAs including a 14-state logically universal NCCA with 
rotation-symmetry [S], always with square neighborhoods. 

In this paper, we show specific necessary conditions for triangular and hexag- 
onal CAs to be number-conserving. Under some symmetry assuptions (rotation 
symmetry for the triangular case and permutation symmetry for the hexagonal 
one) , we show that the local transition function can be decomposed into the sum 
of several flow functions, that is, functions only depending upon two variables. 
These flow functions are later used to design respective simulations between 
hexagonal cellular automata and triangular cellular automata if we assume both 
cellular automata to be permutation-symmetric. 

This paper is organized as follows; section [2] recalls the classical definitions 
that will be used; section [3] exhibits the necessary conditions for triangular and 
hexagonal CAs to be number-conserving. And, finally, in section |4l we present 
number-conserving simulations between hexagonal cellular automata and trian- 
gular cellular automata under the permutation symmetry assumption. 

2 Definitions 

Definition 1. A deterministic two-dimensional radius one cellular automaton 
is a 5-tuple defined by A = [1? ,n,Q, f,q), where Z is the set of all integers, 
n G {3,4,6} is the number of neighbor cells (which implies a corresponding 
neighbor vector set, which is a finite and ordered set of distinct vectors from 
1? : {v^, . . . ,v^), Q is a non-empty finite set of internal states of each cell, 
f : Q" Q is a mapping called the local transition function and q G Q is a 
quiescent state that satisfies /(<?,•••,(?) = q- 

A configuration over Q is a mapping a : ^ Q. The set of all configurations 
over Q is denoted by Coni{Q), i.e., Conf((5) — {a\a : I? Q}. The function 
F : Conf (Q) Conf(Q) is defined as follows and is called the global function oi A 
induced by /: Va e Conf((3), Vlf G Z'^,F{a){l?) = f{a{lf+vo), . . . ,a{lf+vi)). 
From now on, we will denote if by {x, y). 

Let us denote hy Cp the set of finite configurations i.e. which have a fi- 
nite number of non-quiescent states. A cellular automaton A is finite number- 
conserving (FNC) when it satisfies 

VaeCf, {Fia)ix,y)-a{x,y)} ^0. 

And, according to 3J, A is FNC if and only if it is number-conserving. 
Next we define some symmetry conditions of common use (eg. see [5]). 

Definition 2. CA A is rotation-symmetric if its local function f satisfies: 



Vg,Si e Q {I < i < n),f{g, si, ■ • ■ , s„) = f{g, S2, • ■ • , Sn, si), 
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and A is permutation-symmetric if its local function f satisfies: 

yg,S^ e Q {i < 1 < n),V7r e 5„, f{g, Sl, • ' ' , Sn) = ■S7r(l), ■ • • , s^{n)), 

where Sn denotes the permutation group with n elements. 
2.1 Simulation 

Below, we propose the definition of a step by step simulation between two CAs. 
It expresses that if a CA A simulates each step of CA B in r units of time, there 
must exist effective applications between the corresponding configurations [6]: 

Definition 3. Let Conf^i and Conf^ be the two sets of CA configurations of 
(resp.) A and B. We say that A simulates each step of B in time t (and we 

note B < A) if there exists a constant t £ N and two recursive functions k : 
Conis Conf/i and p : CoidA Conf^ such that k o p = Id and for all 
c,c' € Confs, there exists c" € Conf^i such that if cJ ~ Fb{c), cJ' =^ F^{k{c)) 
with p{c") = d , where Fm denotes the global transition of CA M and Fj^ the 
t-th iterate of a global transition of CA M . 

Depending upon the value of t, we say that the simulation is elementary if 
T — 1 and simple if r = 0{1). 

3 Von Neumann neighborhood number-conserving CA 

Durand et al. [2] proved a general necessary and sufficient condition for a NCCA 
with n X m neighbors to be number-conserving. With this condition, any local 
function can be decomposed into the summation of the local function in which 
several arguments are fixed to zero (which is a quiescent state). The drawback 
of this general statement is that it does not explicitely represent neither the 
movement of the numbers nor symmetries. In the sequel, we show novel neces- 
sary conditions for NCCA in different lattices structures, namely triangular and 
hexagonal. The case of the square grid was already considered in [8], where a 
necessary and sufficient condition for a von Neumann neighborhood CA to be 
number-conserving was shown. 

3.1 Triangular number-conserving cellular automata 

Theorem 1. A deterministic two-dimensional rotation- symmetric triangular 
CA A = [1? ,?t,Q,t,q) is number- conserving iff t satisfies 

3ip : ^ Z, \/g. a,b,c ^ Q, 

t{g, a,b,c)^ g + ip{g, a) + ip{g, b) + Lp{g, c) 

fig, a) = -f{a,g)- 



with ip{g, a) = t{g, a, q, q) - t{g, q, q, q) - t{q, g, q, q) + q. 
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Fig. 1. A configuration in the triangular case. 



Proof. Let S{g,a,b,c) = t{a,g,q,q) + t{b,g,q,q) + t{c,g,q,q) + t{g,a,b,c) + 
2t{q, a, q, q) + 2t{q, b, q, q) + 2t{q, c,q,q) — a — b — c — g — 6q. With respect to 
the configuration of Fig[l] only shadowed cells change their states in the next 
step. Then for any g, a, 6, c in Q, 6(g, a, 6, c) = is necessary to preserve number 
conservation. Let's consider the following equation. 

S{g, a, b, c) - S{g, a, q, q) ~ S{g, q, b, q) - 5{g, q, q, c) + 3S{g, q, q,q) ^0 

To satisfy the number-conservation, it is also necessary. Finally, the following 
condition is necessary by expanding the equation. 

t{g, a, b, c) = g+3q+t{g, a, q, q)+t{g, q, b, q)+t{g, q, q, c)-3t{g, q, q, q)-3t{q, g, q, q) 

Let ifig, a) = tig, a, q, q) - t{g, q, q, q) - t{q, g, q, q) + q, then ip{g, a) + ifi{a, g) = 
2q -t{q, g, q, q) - t{q, a, q, q) - t{g, q, q, q) - t{a, q, q, q) + t{g, a, q, q)+t{a, g, q, q) = 
t{g, a, q, q) - t{g, q, q, q) - t{a, q, q, q) = 0. 

We use Durand et al. result ^3] for proving the sufhciency. □ 

Remark 1. The condition also holds in the case of permutation-symmetry. 



3.2 Hexagonal number-conserving cellular automata 

Theorem 2. A deterministic two-dimensional permutation- symmetric hexago- 
nal CA, Ti ~ (Z^, 6, Q, 5, q) is number- conserving iff its local transition function 
5 satisfies: 

3tp:Q'^^Z, yg,a,b,c,d,e, f,g € Q, 

6{g, a, b, c, d,ej) = g + tpig, a) -f ^{g, b) + tp(g, c) + i){g, d) + tp(g, e) -I- ip{g, /), 
with %l){g, x) = S{g, x, q, q, q, q, q) - 6{g, q, q, q, q, q, q) - S{q, g, q, q, q, q, q) + q. 
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Fig. 2. A configuration in the hexagonal case. 

Proof. We show that the condition is necessary. Let us assume that H is FNC; 
then according to Fig. [21 

g + a + b + c + d + f + 12q = S{g, a, b, c, d, e, /)+ 

Ha, b, f, g, q, q, q) + 6{b, a, c, g, q, q, q) + S{c, b, d, g, q, q, q) + 

S{d, c, e, g, q, q, q) + 5{e, d, f, g, q, q, q) + 5{f, a, e, g, q, q, q) + 

S{q, a, q, q, q, q, q) + 5{q, b, q, q, q, q, q) + 6{q, c, q, q, q, q, q)+ (1) 

^{q, d, q, q, q, q, q) + 6{q, e, q, q, q, q, q) + 5{q, a, q, q, q, q, q) + 

S{q, a, b, q, q, q, q) + 6{q, b, c, q, q, q, q) + 6{q, c, d, q, q, q, q) + 

S{q, d, e, q, q, q, q) + S{q, e, /, q, q, q, q) + 6{q, f, a, q, q, q, q) + 

The local function S{g,a,b,c,d,e, f) which satisfies equation ([T]) only depends 
upon terms of four non-quiescent variables. The idea is to decrease this number 
of variables down to two, to finally get the binary ijj function. We first prove 
Lemma [T] which allows to decrease the number of variables in S. 

Lemma 1. In the case of permutation- symmetry, the following equation holds 
for a hexagonal CA, TC = (Z^, 6, Q, S, q): \/g, x,u, z £ Q, 

Hg, x,y,z,q,q,q) ^ g + X + y + z + Uq - d{x, g, q, q, q, q, q) - 6{q, x, y, g, q, q, q) 
-S{y, 9, q, q, q, q, q) ~ S{q, y, z, g, q, q, q) - 5{z, g, q, q, q, q, q) - 6{q, x, z, g, q, q, q) 
-3(5(q, X, q, q, q, q, q) - 36{q, y, q, q, q, q, q) - 3S{q, z, q, q, q, q, q) 

Proof. Cancelling variables 6, d and / by assigning them to the quiescent state 
q in equation ([T]) gives, because Ti. is permutation-symmetric: 

g + a + c + e + 12q = S{g,a,c,e,q, q, q) + d{a, g, q, q, q, q, q) + 5{q, a, c, g, q, q, q) + 
^(c, g, q, q, q, q, q) + S{q, c, e, g, q, q, q) -f 5{e, g, q, q, q, q, q) + S{q, a, e, g, q, q, q) + 
3S{q, a, q, q, q, q, q) + 36{q, c, g, q, q, q, q) + 36{q, e, q, q, q, q, q) 

□ 

Lemma 2. In the case of permutation- symmetry, the following equation holds 
for a hexagonal CA, TC = (Z^, 6, Q, S, q): V.t, y Q, 

6{q, X, y,q,q,q,q) = nq + X + y 

-56{q, X, q, q, q, q, q) - 5S{q, y, q, q, q, q, q) 
-S{x, q, q, q, q, q, q) - (5(y, q, q, q, q, q, q) 
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Proof. Lemma [2] is proved by setting b, c, e and f to q in equation ([T]). □ 

Lemma 3. In the case of permutation- symmetry, the following equation holds 
for a hexagonal CA, Ti. = (Z^, 6, Q, 6, q): \/x, y Cz Q, 

8q + x + y^ 3S{q, x, q, q, q, q, q) + 36{q, y, q, q, q, q, q) 
+6{q, X, y, q, q, q, q) + S{q, y, x, q, q, q, q) 
+5{x, y, q, q, q, q, q) + S{y, x, q, q, q, q, q). 

Proof. Lemma |3] is proved by replacing g,c,d and e by g in equation ([T]). □ 

Lemma 4. In the case of permutation- symmetry, the following equation holds 
for a hexagonal CA, Ti. = (Z^, 6, Q, S, q): Vx, y ^ Q, 

X = -6q + 65{q, X, q, q, q, q, q) + S{x, q, q, q, q, q, q). 
The proof of Lemma |4] is straightforward. 

We now prove Theorem [2j We first make a repeated use of Lemma [1] by sus- 
tracting it from equation ([T]) with suitable variables substitutions and we obtain 
equation ([2]) which only depends upon terms in two non-quiescent variables. 

5(a + 6 + c + d + e + /) + 7g + 174g+ 

S{a, b, q, q, q, q, q) + 5{a, f, q, q, q, q, q) + S{b, a, q, q, q, q, q) + 5{b, c, q, q, q, q, q) + 
S{c, b, q, q, q, q, q) + 5{c, d, q, q, q, q, q) + 5{d, c, q, q, q, q, q) + 5{d, e, q, q, q, q, (?) + 
S{e, d, q, q, q, q, q) + S{e, f, q, q, q, q, q) + 5{f, o, q, q, q, q, q) + S{f, e, q, q, q, q, g) + 
&{g, a, q, q, q, <?, q) + 3{g, b, q, q, q, q, q) + 5{g, c, q, g, q, q, q) + S{g, d, q, q, q, q, q) + 
S{g, e, q, q, q, q, q) + S{g, f, q, q, q, q, q) = 5{g, a, b, c, d, e, /) + 125{g, q, q, q, q, q, g) + 
16((5(g, a, q, q, q, q, q) + S{q, b, q, q, q, q, q) + S{q, c, q, q, q, q, q) + 5{q, d, q, q, q, q, q) + 
S{q, e, q, q, q, q, q) + 5{q, f, q, q, q, g, g)) + 185(g, g, q, q, q, q, g) + 2S{q, a, f, q, q, q, g) + 
26{q, a, g, q, g, q, q) + 2S{q, b, a, q, q, g, g) + 2S{q, b, g, g, g, g, g) + 25{q, c, b, q, g, g, g) + 
2S{q, c, g, q, q, q, q) + 25{q, d, c, g, g, g, g) + 25{q, d, g, q, g, g, g) + 2S{q, e, d, q, g, g, g)-|- 
2S{q, e, g, q, q, q, q) + 2S{q, f, e, g, g, g, g) + 25{q, f, g, q, q, q, q) + 2(5(g, g, a, q, q, q, q) + 
2S{q, g, b, q, q, g, g) + 25{q, g, c, q, q, q, q) + 25{q, g, d, q, q, q, g) + 2S{q, g, e, g, g, g, g) + 
2S{q, g, f, q, q, g, g) + 3(5(g, a, b, g, g, g, g) + 5{q, b, c, q, q, q, q) + S{q, c, d, q, g, g, g)-|- 
S{q, d, e, g, g, g, g) + 5{q, e, f, q, g, g, g) + S{q, f, a, q, q, g, g)) + 7{5{a, g, g, g, g, q, g) + 
S{b, g, g, g, g, g, g) + 3{c, g, g, g, g, g, g) + 3{d, g, g, g, g, g, g) + 5{e, g, g, g, g, g, g) + 
S{f, q, q, q, q, q, q)) + S{q, a, c, g, g, g, g) + 5{q, a, e, q, q, q, g) + 5{q, b, d, q, q, q, g) + 
S{q, b, f, q, q, q, q) + S{q, c, e, q, q, q, q) + S{q, d, f, q, q, q, q). 

(2) 

We observe that we have two kinds of terms in equation 

~ S{q,x,y,q,q,q,q); 
- 5{x,y,q,q,q,q,q). 
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The former will be changed by the repeated use of Lemma [2] and the latter by 
the repeated use of Lemma |3l This yields to equation ([3]) . 

13a + 136 + 13c + 13d + 13e + 13/ + 17g + blOq + 5{g, a, b, c, d, e, /) = 
13d{a, q, q, q, q, q, q) + 13S{b, q, q, q, q, q, q) + 13^(c, q, q, q, q, q, q) + 
13S{d, q, q, q, q, q, q) + 13(5(e, g, q, q, q, q, q) + 13S{f, q, q, q, q, q, q)+ 
5{g, a, g, g, g, g, g) + Sig, b, q, q, q, q, g) + 6{g, c, q, q, q, q, q)+ 
S{g, d, q, q, q, q, q) + S{g, e, q, q, q, q, q) + 6{g, /, q, q, q, q, q)+ 
78d{q, c, q, q, q, q, q) + 78S{q, e, q, q, q, q, q) + 185{q, b, q, q, q, q, q)+ 
78S{q, f, q, q, q, q, q) + 786{q, a, q, q, q, q, q) + 185{q, d, q, q, q, q, g)+ 
^'^S{g,q,q,q,q,q, q). 

The sum of the neighbors 13(a + b + c + d + e + f) is removed by the repeated 
application of Lemma S) Every time we use Lemma SI we also cancel terms of 
the form 785{q, x, q, q, q, q, q) + 135(a:, q, q, q, q, q, q) in the rhs of equation [3l All 
remaining terms in the rhs of equation [3] are like S{g, x, q, q, q, q, q). 

We use Durand et al. result [3| for proving the sufficiency. This proves The- 
orem [2l □ 

All above computations were made using computer algebra systems and the 
spreadsheets can be obtained from the authors. 

4 Simulations between hexagonal and triangular NCCAs 

In this section, we propose effective mutual simulations between triangular and 
hexagonal permutation-symmetric NCCAs. 

4.1 Elementary simulation of a triangular NCCA by a hexagonal 
NCCA 

Proposition 1. For any triangular permutation- symmetric NCCA T , there is 

1 

a hexagonal permutation- symmetric NCCA 7i such that T ^ Ti. and whose tran- 
sition function is surjective. 

Proof. Let T = {1?, 3, Qt, tr, q) be a triangular permutation-symmetric NCCA 
with flow function Lp. 

We construct TL = (7? ,6,QH,tnT<l') a permutation-symmetric hexagonal 
NCCA by designing its flow function 

Let q' ^ Qt and Qn = Qt U {g'}; the flow fimction corresponding to t-^ 
contains the following values: 

For each x,y G Qt, assign 'ip{x, y) = 1^9(2;, y), 
For each x G Qr, assign , ipix, q) = 0. 

Given an initial configuration of T (see Fig. [S]), let the initial configuration 
of H be as depicted in Fig. El From the process of Ti, it is clear that the cells 
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with value q' don't have any effect and the other non-quiescent cells simulate T. 
Applying the algorithm in 5 , it is possible to add extra-states produced by tu 
to Qt-i and make the local function tji to be surjective. □ 



\Cb.3 



\C2,3 



\Cb.i 



,0 \ / »-'l,0 ^ 



\Cb.i 



\Co,- 



/Cb.-: 



/Ca,-2\ /Ci,-: 



Fig. 3. An initial configuration of triangular CA T. 




4.2 Simple simulation of a hexagonal NCCA by a triangular NCCA 
Proposition 2. For any hexagonal permutation- symmetric NCCA Ti., there is 

2 

a triangular permutation- symmetric NCCA T such that Ti. ^ T and whose tran- 
sition function is surjective. 

Proof. Let TL = (Z^, 6, Q^^, g) be a hexagonal permutation-symmetric 
NCCA. We assume that Q = {sq, si, • • ■ , Sm-i} and its flow function is ■(/; and 
Ti. has an initial configuration as pictured in Fig. O 

Without loss of generality, we can assume < min{si) < max{sj) < M for 
a constant M, because non-positive numbers in Q-j-c can be erased by adding a 
constant value to every state numbers and by changing the arguments of the 
rules. 
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Fig. 5. An initial configuration of hexagonal CA TL. 



We construct T = (Z^ 3, Qr, ^r, 0) a permutation-symmetric triangular 
NCCA by designing its flow function Lp. 

First we assign numbers corresponding to each state in Qu- Let pi = 
2rigMl+,(^^0,---,m-l). 

The flow function ip of tr contains the following values: 

For each Si S Qu, assign 

</j(0, Si) = Pi. 

For each combination of Sj, Sk G Qh, 

ip{-3p^,p^ +pj +pk) ^Pi + ip{si,Sj) +i>{st,Sk)- 

All other values of cp are 0. The local function tq- can also be extended to be 
surjective. 

The initial conflguration of T is chosen as in Fig. [6l 

We briefly explain how the rules work. In the first step, each nonzero state 'x' 
in Fig. [S] moves three values Px to the neighboring zero cells; after this, the zero 
cells contain the value Px+Py+Pz related to the three neighboring cells ^x\^y\'z\ 
In the second step, the cell on which we focus (was 'x') knows the values of two 
neighboring cells of Ti. by the neighboring cell which value is Px + Py + Pz and 
can move the values of tpisx, Sy) + ipisx, Sz)- □ 

5 Conclusion 

In this paper, we have designed flow functions for number-conserving triangular 
and hexagonal cellular automata under the permutation-symmetry condition. 
This was also generalized to the rotation-symmetry for triangular NCCA. A 
simulation between triangular and hexagonal NCCA, and conversely, were also 
proposed. 
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Fig. 6. An initial configuration of triangular CA T. 



This work can be extended in several ways. First, we'd like to know if the 
flow function we proposed in the permutation-symmetry hexagonal case also 
holds in the rotation-symmetric case. Then, we aim to go on in the number- 
conserving simulation of different neighborhoods like simulating a von Neumann 
neighborhood in a square lattice by a triangular NCCA. 

Finally, it might be possible to generalize these results to the case of number- 
conserving cellular automata on Cayley graphs and to find flow functions in 
several cases and number-conserving simulations as well. 

Acknowledgments. This work was done while K. Imai was visiting I3S labo- 
ratory thanks to a grant from the french CNRS institution. 
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